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Abstract 

We performed first principles calculations on work functions of single wall carbon nanotubes 
(SWNTs), which can be divided into two classes according to tube diameter {D). For class / tubes 
{D > Inm), work functions lie within a narrow distribution (~ O.leF) and show no significant 
chirality or diameter dependence. For class // tubes {D < Inm), work functions show substantial 
changes, with armchair tubes decrease monotonously with diameter, while zigzag tubes show op- 
posite trend. Surface dipoles and hybridization effects are shown to be responsible for the observed 
work function change. 
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Carbon nanotubes (CNTs)^ have attracted considerable attention as a promising can- 
didate material for future nano-scale electronics. It boasts reliable structural strength, ma- 
ture synthesis methods, and unique electronic structure that is closely related to its chirality. 
Nanotube with indices (n,m) can be either metallic or semiconducting depending on whether 
n — m is a multiple of three. Besides bandgap, absolute potential of the fermi level, namely 



unc- 



the work function, is another important physical quantity of interest. Nanotube work 
tion is a critical quantity in understanding field emission properties from SWNTsj^. It 
dictates the direction of charge transfer and influences the selectivitv of molecules as chem- 
ical sensors It affects band lineup at the CNT/metal contact p| and has considerable 
impact on device performance j^. From device integration perspective, it is desirable to 
have all the nanotubes acting as channel materials with the same work function. All these 
indicate the importance of a thorough understanding of nanotube work functions. 

However, there have only been few theoretical studies on the work functions of SWNTs. 
Simple TT-orbital tight binding, which serves as a good starting point for understanding 



nanotube electronic structure, would predict no work function change at all. Zhao et al 
studied work functions of CNTs based on the molecular orbital calculations, indicating 
higher work function of semiconducting nanotubes and considerable diameter dependence 
over all diameter range. Some other ah inito calculations ^, 0, [3| have investigated cur- 
vature effects and their effects on the band gap. Less is clear about their implications on 
work functions. On the experimental side, it has recently become possible to measure work 
functions of individual SWNTs However, there is so far no consensus among dif 



ferent experiments. Resonant confocal micro Raman spectroscopyjlll] reveals considerable 
work function difference between metallic and semiconducting nanotubes while photoemis- 
sion electron microscopy study^^ shows work functions of various nanotubes is generally 
insensitive to either diameter or chirality. In this paper, we employed first principles plane- 
wave pseudo potential calculations and systematically studied work functions of (n,0) and 
{n,n) nanotubes. The nanotubes we studied range from ultra small tubes ~ OAnm to con- 
ventional size tubes ~ 2.5nm usually grown by chemical vapor deposition (CVD). Work 
function study of ultra small tubes is motivated by recent experimental fabrication and ob- 
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16|. From a systematic 



servation of nanotubes with diameters down to ~ OAnm 
analysis of the calculation results, work function behavior of nanotubes can be divided into 
two classes. Work functions of class / tubes fall in a narrow distribution ~ O.leV. For 
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class // tubes, work functions increase or decrease significantly depending on whether it is 
armchair or zigzag tube. 

Workfunction is defined as WF = — Ef, where is the vacuum level and Ef is the 
Fermi level of the system. In the case of semiconducting nanotubes, we place Ef at the 
midgap, which is consistent with experimental measurements jlll. 1^ . The work function 
calculation is done using 

the work function, the nanotube is put into a supercell of 35A x 35A in the cross sectional 
direction. This corresponds to a nearest inter-tube distance of > 10 A for the largest tube 
under study. Kohn-Sham single-electron wavefunctions are expanded by 921984 plane waves 
with an energy cut off of 21.06 Ryd. 30 k-points are used along the one dimensional Brillouin 
zone which are tested to give good convergence. All the atoms in the unit cell are fully relaxed 
until the force on each atom is less than 0.02e\^/A. For class / tubes, there is little difference 
between the optimized geometry and rolling up of a graphene sheet. For class // tubes, the 
main effect of relaxation is the expansion of the tube diameter by a few percent Q]. The 
vacuum level (p is determined from the average potential in the vacuum region where it 
approaches a constant. Local density approximation (LDA) are used throughout for the 
exchange-correlation function. It is known that LDA utilized with Kohn-Sham scheme 
tends to underestimate the bandgap, which in our case, results in an over estimate of the 
work function for semiconducting tubes. For an estimate of the error that may arise from 
band gap underestimate for semiconducting nanotubes in class /, we compare our calculated 
LDA bandgap with Tight binding(TB) and GGA calculations. Table I lists the bandgap 
for selected CNTs. It can been seen that for reasonably large diameter tubes, the band 
gaps derived from different theoretical models differ by less than O.leV. This translates to 
an error ~0.05eV in work function and will not affect our conclusions. Also, (r;,,0) tubes 
with n < 6 become metallic due to strong curvature effects and we expect LDA to give a 
reasonably accurate description of their electronic structures. The largest possible error in 
work function estimate then occurs for (7,0) and (8,0) semiconducting tubes, which belong 
to class // tubes. One way to correct the error is to use Hedin's GW method |2o| which has 



2l| . GW calculations 



been shown to give accurate results for a wide range of materials 
confirm the LDA result that (6,0) CNT is metallic and (7,0) is semiconducting (Table I). 
Thus corrections from GW method will only shift work function values for (7,0) and (8,0) 
tubes without altering the overall trend of work function change. 
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FigH] summarizes the work functions of (n,0) and (n,n) nanotubes with diameters ranging 
from 0.32nm to 2.5nm, plotted against the inverse tube diameter. The dotted horizontal line 
in the graph indicates the calculated work function for a graphene sheet of A.66eV. For class 
/ nanotubes, work functions of both semiconducting and metallic tubes fall within a very 
narrow distribution around 4.66eV, just as would be predicted from 7r-orbital tight binding 
model. More detailed analysis shows the work function decreases slightly with diameter. If 
work functions of class / tubes are fitted to a curve of the form y = a/D + 4.Q6eV, then 
a = —0.2 ± 0.08 eV ■ A. The slope of the decrease is so small that for practical purposes, 
work functions of class / tubes can be well approximated by the graphene work function. 
We can conclude that for class / nanotubes, there is no significant work function difference 
within LDA framework. Since the tubes used in both experiment s ll ll . Il2 1 are larger than 
Inm in diameter, our calculation is more consistent with Suzuki's finding that there is 
no structural dependence on the work functions of CNTs. The discrepancy between our pre- 
diction and that of Zhao'sj^] regarding work functions of semiconducting nanotubes largely 
lies in the difference of work function definition, where they approximated the work function 
as WF = —Ehomo for semiconducting tubes. Thus the larger work functions of semi- 
conducting tubes and inverse diameter dependence in Zhao's study is primarily due to the 
changes in band gap. Comparing with their results WF = —Ej for metallic nanotubes, we 
predict an even weaker dependence on diameter for class / tubes by including contribution 
from vacuum level. Besides these differences, our calculation for class / tubes is generally 
consistent with the results from molecular orbital calculations. However, significant changes 
in work functions can be observed for class // tubes. As can be seen from Fig. ^ work 
functions of class II nanotubes split distinctively into two groups. Work functions of {n,n) 
tubes monotonously decreases with diameter, while work functions of (n,0) tubes show a 
dramatic increase. The magnitude of the work function change is comparable to nanotube 
bandgap and thus plays an important role in determining the band lineup [3]. Its conse- 
quences can be experimentally observed via transport measurements!^, 24 1. For the smallest 
tube we studied, (4,0) tube has an ultra high work function of 5.95eV^. 

Work function increase (decrease) can be induced by either an enhanced (reduced) sur- 
face dipole, or a lowering(rising) of its intrinsic bulk Fermi energy. To gain a quantitative 
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understanding from these two contributions, we rewrite work function as 

WF = <P-Ef=[cP- Vref] - [Ef - Vref] (l) 

where Vref represents average electrostatic potential at carbon atomic core. In this expres- 
sion, the first term is the potential with respect to the vacuum level (differ by a minus sign), 
whose change reflects a change in surface dipole contribution as we outline below. The sec- 
ond term represents the relative position of the Fermi energy above the average potential, 
which is related to the material's intrinsic bulk electronic structure. As shown in FigI2l 
the average potential at the carbon atom monotonously decreases with decreasing diame- 
ter. By comparing with the average potential in a graphene sheet of corresponding periodic 
boundary condition, it can be concluded that the lowering of the potential is not due to 
quantization along the circumference in a nanotube, but rather due to curvature effect. For 
large diameter tubes, the potential inside and outside the tube are almost equal, indicating 
symmetric charge distribution. As the tube becomes smaller, the potential inside the tube 
falls below the vacuum level outside, indicating asymmetric charge distribution, with more 
charge distributed outside. More asymmetric charge distribution results in an enhanced 
surface dipole and lowers Vref with respect to vacuum level. The potential lowering due 
to surface dipole is insignificant for class / tubes, but it can be up to 0.5eV for class II 
tubes. The first term in Eqn. (0) alone, however, does not dictate the trend of work function 
change. Accompanying the lowering of the average potential is the tighter confinement of 
wave function along the circumferential direction, and the electron fermi level is generally 
expected to rise with respect to Vref, which counteracts the effect of work function increase 
due to potential lowering. In fact, for class / nanotubes, the second term increases faster 
than the first term, resulting in a slightly decreasing work function as the tube diameter 
decreases. 

The above picture holds true for class / tubes. However, the second term Ef — Vref 
decreases rapidly for zigzag nanotubes smaller than (10,0). Combined with the potential 
lowering from the first term, this makes zigzag tubes' work functions increase dramatically. 
This is exactly where a singly degenerate state in the conduction band moves down below 
the conduction band bottom due to a* — n* hybridization j^. In Fig. El , we show band 
structures of (10,0), (7,0) and (4,0) CNTs respectively. In (10,0) CNT, the conduction band 
bottom is doubly degenerate as is predicted by simple 7r-orbital tight binding model. The 
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singly degenerate state (indicated by a solid line in the figure), which originates from the 
band that crosses F in the Brillouin zone lies approximately leV above the conduction band 
bottom. Upon increasing curvature, this singly degenerate bands is significantly lowered. 
In (7,0) CNT, it is below the n* conduction band bottom and lies within the bandgap, 
whereas for (4,0), (5,0) and (6,0) CNTs, it is even lowered below the tt valence band top, 
transforming them into metallic tubes. The lowering of this singly degenerate band due to 
curvature effects has been observed in many previous theoretical studies j^, 0, however, 
its implication on work function change has not been examined. Most of the earlier works 
have focused on its effect on band gap change, and the energy levels are studied with respect 
to the fermi level, which is usually set to zero. 

In the following, we show that the lowering of this singly degenerate band has important 
effects on work function. Fig. shows the graphene band structure along F — M line folded 
at the middle of Brillouin zone and energy levels of the corresponding singly degenerate 
bands at F point in a (n,0) nanotube. All the energy levels are drawn in absolute scale with 
respect to vacuum (which is set to be zero). The energies of these bands at F point show a 
downward shift compared to that of graphene, indicating the lowering of the potential inside 
the tube (Fig. The top of the singly degenerate tt band is lowered ~ 1.0e\^ in going from 
graphene to (4,0) CNT, and a band is lowered by ~ 0.5eV, roughly in accordance with the 
lowering in K-e/- The effect of a* — it* hybridization induced by curvature is to bring down 
the singly degenerate vr* band with respect to the already lowered a and tt bands. Before it 
moves below the doubly degenerate conduction band bottom, the work function is largely 
determined by the doubly degenerate valence band top and conduction band bottom. Once 
this singly degenerate states is down shifted below the conduction band bottom, as is the 
case for (7,0)and (8,0) CNT, there is an accompanying increase in work function due to 
the lowering of the conduction band edge. For zigzag tubes smaller than (6,0), this singly 
degenerate state crosses the valence band, pinning the fermi level below the valence band 
top. The concerted effects of potential lowering and moving down of the singly degenerate 
state makes the (4,0) CNT a material with extremely high work function. Fig. |3Jd shows 
the same plot for {n,n) tubes, a band top is also lowered by similar amount as K-e/? but due 



26|, there is no appreciable lowering of the 



to weaker a* — tt* interaction in armchair tubes 
singly degenerate ir* state. In fact, more pronounced effect is the upshift of the bonding vr 
orbital due to enhanced a — n interaction ^ which results in an decrease in work function. 
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In conclusion, we have studied work functions of intrinsic SWNTs with experimentally 
observable diameter ranging from ~ OAnm to ~ 2.5nm. Work functions of class / tubes 
show little dependence on diameter or chirality. For class // tubes, work functions show 
significant differences for {n,0) and {n,n) tubes, a* — n* interaction which pushes a singly 
degenerate state to cross the valence band is the main reason for the dramatic increase of 
work function observed in (n,0) tubes. 
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Chirality TB* GW[22] GGA[25] Our 


Data(LDA) 


(6,0) 0.14 metal 


metal 


metal 
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0.19 
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(10,0) 0.85 / 


0.88 
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0.61 


0.57 



*tight binding parametrization is similar to that in Ref| 
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FIG. 1: Work functions of {n,n) and (n,0) nanotubes of different diameter. For class / tubes, 
nanotube work functions are very close to the work function of a graphene sheet ~ 4:.66eV. For 
class II tubes, work functions of (n,0) tubes increase dramatically, while those of tubes show 
decrease in work function. (6,0), (5,0) and (4,0) tubes are metallic due to a* — vr* hybridization. 
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FIG. 2: As tube diameter decreases, there is a lowering of the average potential at carbon atomic 
core with respect to vacuum which is due to asymmetric charge distribution inside and outside the 
tube. This effect is more pronounced for class // tubes where a slight change in diameter leads to 
a large reduction in potential energy. The graphene sheet with corresponding boundary condition 
does not show such reduction. 
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FIG. 3: Band structures of (a) (10,0) (b) (7,0) (c) (4,0) CNTs respectively, here a = 3ao, where 
gq = 1.44A is the lattice constant of graphene. The singly degenerate state (solid line) lies above the 
conduction band bottom (a), in the bandgap (b) and below the valence band top (c) respectively. 
Accompanying the lowering down of the singly degenerate state is the lowering of the Fermi level 
with respect to the vacuum. 




FIG. 4: Evolution of the singly degenerate energy levels at F point for (a) (n,0) and (b) (n,n) 
CNTs. On the right hand side of both graphs is the graphene band structure folded at the middle 
of Brillouin zone, where a = Sgq and b = \/3flo- On the left hand side of both graphs is the 
absolute energy values of singly degenerate states in the nanotube, which is derived from the band 
that crosses F point in graphene in the zone folding model. 
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